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. Abstract 

We discuss the general-relativistic contributions to occur in the electromagnetic 
properties of a superconductor with a heat flow. The appearance of general-relativistic 
. contribution to the magnetic flux through a superconducting thermoelectric bimetallic 
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observable in the superconducting state in the Earth's gravitational field are considered. 
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circuit is shown. A response of the Josephson junctions to a heat flow is investigated 
in the general-relativistic framework. Some gravitothermoelectric effects which are 
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1 Introduction 



Effects of fields of gravity and inertia on a superconductor have been investigated by a 
number of authors starting DeWitt [1] and Papini [2]. More recently, a general relativistic 
treatment for electromagnetic effects in normal conductors with a heat flow and superconduc- 
tors without gradient of temperature has been given. In particular, several superconducting 
devices that can, in principle, detect a gravitational field have been presented by Anandan 
[3]. Further, one of them has been already realized in the experiment [4] which tested the 
equivalence principle for Cooper pairs. 

From our point of view a study of thermoelectromagnetic relativistic gravitational effects 
in superconductors with nonzero gradient of temperature is of fundamental interest due to 
the following two reasons. 

Thermoelectric effects do not vanish, in principle, in inhomogeneous and anisotropic su- 
perconductors and nowdays attract an increasing interest [5] especially due to the discovery 
of high temperature superconducting materials which combine both anisotropy and inhomo- 
geneity properties and which, in this connection, are favourable for measurement of thermal 
effects. Since superconductors provide extremely sensitive and accurate measurements, there 
is hope that the very weak general relativistic contributions to the thermoelectromagnetic 
effects due to the inhomogeneity arising from the Earth's gravitational field might be de- 
tectable. 

On the other hand, according to the recent theoretical models, core of neutron stars 
forms matter in a superconducting state with thermal distribution and heat flow (see, for 
example, [6]). It suggests that thermoelectromagnetic effects might explain of origin and 
evolution of magnetic field inside the core of supermassive pulsars for which the dimensionless 
general relativistic compactness parameter a/R a reaches about 0.5 [7] (a and R s are the 
gravitational radius and radius of star, respectively). In addition thermoelectromagnetic 
effects in Josephson contacts could be responsible for electromagnetic radiation arising from 
possible superconductor-normal metal-superconductor layers in the intermediate boundary 
between conducting crust and superconducting core inside neutron star. 
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2 General Relativistic Correction to the Magnetic Flux 
Induced by a Heat Flow in Superconducting Ther- 
mocouple 

Consider a superconductor with a heat flow in the external stationary gravitational field. 
According to the two-component model two currents flow in this superconductor: the su- 
perconducting current of density jr s \ a and the normal current of density j( n )a- The normal 
current is carried by 'normal' electrons (exitations) and it does not differ essentially from 
the current in the normal state of a metal [8] 

F aP U P = jj {n )a + R H (F ua + u a u° F ua )f (n) + A' 1/2 V a jl- PA' 1/2 V« T - bf {n) A a p, (1) 

where f = A 1/2 T,/i = A 1 / 2 //, jj, is the chemical potential per unit charge, T is the tem- 
perature, A is the conductivity, A = — £ Q is a timelike Killing vector being parallel 
to the four velocity of the conductor u a , f3 is normal differential thermopower coefficient, 
F a/ 3 = Ap^ — A a £ is the electromagnetic field tensor, A a is the vector potential, Rh is 
the Hall constant, Ap a = U{ a ,p\ + u\pw a j is the relativistic rate of rotation, w a = u^pu 13 is 
the absolute acceleration and b is the parameter for the conductor, \7« and [• • •] denote the 
spatial part of covariant derivative and antisymmetrization. 

Thus the electric current flowing in the conductor in the general case arises for the 
following reasons (with general relativistic corrections and contributions) - (a) electric field, 
(b) Hall effect, (c) nonequilibrium effects and (d) Coriolis (gravitomagnetic) force effects 
described by the last term on the right hand side of equation (1). 

If the wavefunction of the Cooper pairs is tp — n 1 J 2 e l ' & then four-vector of supercurrent 
density is [3] 

ieh 2e 2e 

j{s)a = — {^*(d a - ij-AaW - i>(d a + i—A a )^*}, (2) 



and satisfies to equation 



2n s e r . , 2n s e 



j(s)a = —{-hd a d + 2e/cA a } = -^P a (3) 
m s m s 

where n l J 2 represents the density of Cooper pairs, m s and P a = p a + ^f^4« are the mass 

and the generalized momentum of the Cooper pair and $ is the phase of superconducting 

wavefunction. 
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Since any supercurrent flows only on the surface within the penetration depth, in the 
interior of the superconductor supercurrent (3) is parallel to the four-velocity u a i.e. P a = 
— ^\iu a . It follows from it that interior of the superconductor 

and after multiplying it by u a one can get that [3] 

Ep - A~ 1/2 dpjl = (4) 

everywhere inside the superconductor in the steady state. 

Formula (4) can be also derived from the general relativistic London equations [9] 

2n s e 2 

3(»W<*] ~~ i(s)[a( lnn s),/3] H —F a p - 2cen s (A a p + u [a wp]) = 0, (5) 

which have been obtained by requiring, that inside a superconducting medium, the motion 
of the Cooper pairs is free of resistance. 

Suppose that two ends of a bulk piece of superconductor are at different temperatures, 
T\ and T 2 . The temperature gradient will produce a force on normal exitations of the 
superconductor initiating a current of the normal exitations 

j (n)a = \- l/2 \(3d a f - \R H {F va + u a u a F ua )f {n) + \bf (n) A a p (6) 

as a consequence of Ohm's law (1) under the conditions corresponding to equation (4). 

It therefore follows that below T c under steady-state conditions when equation (4) should 
be obeyed, if d u T ^ 0, the density of the normal current should be finite because of 
equation (6). However, if the circuit is open, then the total current is zero and in the simplest 
case the density of the total current also vanishes: j a = j( s ) a + ]( n ) a = 0, i.e. the normal 
current density j( n ) is cancelled locally by a counterflow of supercurrent density j( s ) 

j(s)a = -J(n )a = -A^X^T + \R H {F va + U^if F va )f (n ) ~ ^f {n) A a p. (7) 

Suppose the superconductor is embedded in the Schwarzschild space-time 

ds 2 = -(1 - a/r)(dxo) 2 + (1 - a/ry'dr 2 + r 2 {dQ 2 + sin 2 6d V 2 ), (8) 

where timelike Killing vector £ a can be chosen so that A = 1 — a/r. If curvature effects are 
negligible, than the apparatus may be regarded as having an acceleration, g, relative to a 
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Figure 1: Totally superconducting circuit made of two metals I and II in the Earth's vertical 
gravitational field. The magnetic field $t through the ring has gravitational contribution 
proportional to gH/c 2 . 

local inertial frame, and thus, A = (1 + 2gH/c 2 ) 2 , where H is the height above some fixed 
point. 

Because of the cancellation of the thermoelectric current in the superconductor, schemes 
for measuring the thermoelectric effects, based on inhomogeneous or anisotropic supercon- 
ductor configurations, become necessary. A review of various experiments, as well as a 
bimetallic superconducting ring has been given e.g. by Van Harlingen [10]. 

We shall now consider two bulk specimens of two dissimilar superconductors, Sj and Su, 
which are brought into contact in such a way that together they form a closed ring, as in 
Figure 1. Suppose that the temperatures of the upper and lower contacts are kept at the 
different temperatures Ti and T 2 respectively. The presence of an open gap in a massive 
circuit when its thickness d (for example, the diameter of the wire forming the circuit) is 
much bigger than the depth of penetration 5 of the field allows us to calculate the magnetic 
flux $ 6 across the gap without solving the problem completely. As it was shown above in 
the bulk of a superconductor j a = + ]( n ) a = and so 

j {n> = -K-^\(3d v f = - 3{s)u = ^[hd v $ - -A v \. (9) 

m s c 

For the sample embedded in the Schwarzschild space-time last two terms in formula (7) 
disappear since in this case magnetic field does not penetrate inside the superconductor and 
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all components the relativistic rate of rotation are zero. 

Integrating the equation (9) along the contour C, which is in the bulk of the superconduc- 
tor and using that § A a dx a — \ J F a pdS al3 = $t and § d v -&dx v = 2irn, where n — 0, 1, 2, ... 
yields directly 

$ 6 = n$ + [ A-^X^f, (10) 

where m s /e 2 n s = A = Atc5 2 /c 2 and $ — nhc/e = 2 x 1CT 7 Gauss ■ cm 2 is quantum of 
the magnetic flux. The current I s which leads to the appearance of a flux <3> fe flows on the 
internal surface of the circuit in a layer of thickness of the order of S. 
If A = A(0) at then from (10), the magnetic flux through contour is 

$t = n$ + A(O)- 1 / 2 ^^ h\d u fdx u = 

Til C 

n$o + - (0A)«)PU1 + gHjc 2 ) - T 2 (l + gH 2 /c% (11) 

where f3i and f3u are the values of (3 for the two metals, H\ and H2 are the heights of the 
junctions above the Earth's surface. This is the general relativistic generalization of the 
thermoelectric effects in the inhomogeneous (bimetallic) superconductor. 

If, for the sake of simplicity, we assume that (f3\)n ^> (f3\)i we then find from equation 
(11) that 

When apparatus is horisontal, the magnetic flux in the circuit would remain flat space- 
time value and last term in formula (12) will disappear. In this case equation (12) describes 
unquantized thermoelectric flux $r in the presence of a temperature gradient. For typical 
parameters of the superconductor the thermal-current-related flux (first term in (12)) is 
expected to be of order of 10~ 2 <3> . In actual fact, however, much - indeed orders of magnitude 
- stronger fluxes of hundreds of $ have been measured experimentally [11]. 

However if the apparatus is brought into a vertical plane, then the magnetic flux will 
be changed according to formula for The height between two junctions is changed by 
H — Hi — if 2 and the magnetic flux undergoes a fractional change gH/c 2 . Then the general 
relativistic contribution in $ when H = 10cm, in the Earth's gravitational field, will be 
proportional to small dimensionless parameter 10~ 17 . 
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There is no doubt the flux production mechanism discussed in this section may be sig- 
nificant and relevant to the problem of origin and evolution of magnetic fields in isolated 
neutron stars since their radius can be only 1.4 — 3.5 times larger than gravitational radius a 
and their substance may be superconducting of //-type or superfluid at high densities (see 
e.g. [12]). 

3 Thermoelectric Effects in SNS Junctions in the Ex- 
ternal Gravitational Field 

In this section we shall first consider the behaviour of an superconductor-normal metal- 
superconductor (SNS) junction when its S electrodes have different temperatures or, in 
other words, when there is a heat flow through the junction placed in a gravitational field. 
After that, we shall discuss the thermoelectric effects arising when, in addition to carrying 
a heat flow, the junction is placed in both magnetic and gravitational fields. 

Let us suppose that there is a temperature difference AT between superconducting elec- 
trodes of the SNS Josephson junction placed in the external gravitational field (8) (see 
Figure 2). Due to the Kirchoff's first law 

3{s)a = 3{n)a- (13) 

According to Ohm's law (1) a normal component of a current in the junction is 

j (n)a = \E a - \A~ 1/2 d a fl + (3\A~ 1/2 d a f. (14) 

Density of the superconducting current flowing through the junction is related to the 
phase difference = A$ across the junction by 

3(s)a = 3(c)a sin (15) 

where j( c ) Q is the critical value of electric current density. 
Using formulae (13)-(15) and Josephson equation [3] 

_ fi ^ = !£,_2£.4^" ! (16) 
Or c c 

one can obtain 



3(c) 



sin i 



AA- 1/2 {-^(^^A^)+^ a f} 



2edr 



Figure 2: Equivalent circuit for SNS junction with heat flow: j T = \(3A l l 2 gradT is the 
current due to the heat flow across the normal layer. 

which after integration on n a dS will give 

A(O) 1 /^/ (c)Sin0= -|~A 1/2 + /3Af. (17) 

Here R — J is the resistance of the normal layer with length dl, J( c ) = J j( c ) a n a dS is 
electric current, n a is normal vector to the cross section of wire dS, d^jdr = <fi :a u a . 

Therefore if thermoelectric current exceeds the critical current of the Josephson junction, 
then as a consequence of (15) and (17), an alternating current (ac) of frequency 

u = ^-PA-^AT (18) 

is produced and the junction emits radiation with the frequency uj measured by an observer 
at rest with respect to the junction. Formula (18) is the general relativistic generalization of 
thermoelectric ac Josephson effect according to which a temperature difference AT across 
the SNS junction results in electromagnetic radiation, which has been predicted in [13] and 
experimentally confirmed [14]. 

Now we would like to show with help of simple example that general relativistic ther- 
moelectric effects may have astrophysical applications. Suppose that SNS junctions can 
be realized in the intermediate boundary between conducting crust and superconducting 
core inside neutron star. Then as a consequence of equation (18), that is, of the existence 
of thermal analog of the ac Josephson effect, we can predict a new mechanism for electro- 
magnetic radiation production from pulsars: a steady heat flow through SNS junctions in 



intermediate boundary gives rise to Josephson radiation with frequency u. Let us roughly 
estimate the frequency of expecting radiation. If thermoelectric power of conducting crust 
P(T ~ 10 5 °K) ~ 10~ 2 V/K, then cu « 3 x 10 13 A~ 1/2 (T 2 - T X )K, and at a temperature 
difference lCT 7o K we have w~3x 10 5 iJ,2. 

According to (18) the frequency of the junction depends from altitude in gravitational 
field and hence frequences iO\ and w 2 of the junction at Hi and H2 (H = H2 — Hi ; Hi < H2) 
are connected through uj{H 2 ) = u(Hi)(l — gH/c 2 ). If (5 m lO" 5 ^ • K^ 1 and H = 10m, then 
the general relativistic red-shift for the frequency is 

Ugr{Hz) = ^-^-[3 AT « 3 x 10~ 3 AT(K). (19) 

This quantity is weak since it is experimentally very difficult to obtain a great temperature 
difference at the superconducting system. If the radiations are detected by one device, there 
would be no frequency difference between them due to the gravitational redshift effect for 
frequency which will compensate correction (19). Hence this null experiment would confirm 
that the temperature difference, AT 2 and ATi, at the upper and low altitudes are connected 
through 

ATi = AT 2 (1 - gH/c 2 ) (20) 

and varying with height. 

Consider two dissimilar SNS junctions, separated by a height H and connected in par- 
allel by superconducting wires to a common heater source. Suppose that flu and fli are 
the thermoelectric power of the upper and lower junctions, respectively. Now we integrate 
equation (3) over the contour which passes through the interior of superconducting ring with 
two Josephson contacts. Then 

he 

$ 6 = n$ + —(07/-0/), 

where 077 and 07 are the contributions due to the phase discontinuities at the Josephson 
junctions. 

The rate of change of magnetic flux is related to Au = uu — ujj by 

d&b <in^ frc^ 
dr dr 2e 
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so if Auj is not equal to zero, a magnetic flux A<3> b 7^ will be induced. As long as A$ b < $ , 
n will remain constant and A$b will increase linearly with time until A$& = $ 0! then the 
order of the step n will change as flux quantum enters the loop. Thus due to the effect of 
gradient of temperature on the junctions with unequal thermoelectric powers is equivalent 
to having a time dependent flux, as given by the last term on the right hand side of the last 
equation. Then using equations (18) and (20) one can derive that the change in magnetic 
field inside the circuit during the time interval [0, r] is 

A$ fe = c HaT^j - AT 2 (3 n (l - 9 -^)}dr + An$ = 
jo cr 

c [ T ATM - /3 n )dT + An$o- (21) 
Jo 

Thus this particular loop is sensitive to the frequency and in this connection to the 
thermoelectric power difference between the junctions. The independence of magnetic flux 
(21) from gravitational field g confirms the validity of formula (18) for the temperature, 
that is in a gravitational field, during thermal equilibrium, the related quantity T (rather 
T) is constant along the sample. In addition a new experiment, in which the thermoelectric 
response creates flux (21) changing with time will give one more possibility of measurement 
of thermoelectric effects in superconductors. It is more or less important, not only from 
general-relativistic point of view, but also for the new proposals [15] for confirmation of 
some aspects of thermoelectric transport theory. 

When the current exceeds the critical value potential difference V = (3 AT appears across 
the junction due to the thermoelectric effects. Since the thermoelectric power of the junction 
77 differs from that at the junction J, the potential differences across the first and second 
junctions, Vn and V}, respectively, will differ so that AV = Vn — Vi — (Pn — f3i)AT. 
The basic technique for the detection of extremely small voltage differences between two 
Josephson junctions by monitoring of magnetic flux change was firstly developed by Clarke 
[16]. 

In the absence of any additional effects on the Cooper pairs, one would thus expect the 
net EMF in the loop containing the junctions to be (Pn — Pi) AT ~ 10~ n V for the typical 
values of parameters (Pn — Pi) ~ 10~ 6 V/K and AT ~ lO -50 ^. For the loop of inductance 
L the evolution of magnetic field is approximately governed by law ^ = —L^. In this 
connection a nonvanishing value for AV would lead, according to (18) and (21), to a time 
varying current (from zero to the critical maximum value in the range of one number of 
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the step n): j^h = — j^AV, which will induce the above discussed flux A$ b — cj AVdr 
through the loop in the linear regime. 

For measurement of gravitational contribution we can propose the following change. Con- 
nect the Josephson junctions to two independent heaters which have different temperatures, 
such that AT2 = AATi and A is constant. In this case the flux (21) takes form 

A$ 6 = An$ + c f AT^pj - Ap n )dr + Ac f AT^n^-dr. (22) 
Jo Jo c 2 

If the apparatus is horisontal then H = and therefore the last term in the magnetic flux 

change should disappear. When the coil is vertical the last term will increase the rate of 

change of flux according to (22). By detecting this contribution, the gravitational corrections 

to thermoelectric effects can be measured. Taking Af3 ~ 10 _7 cm 1 / 2 ■ g 1 ^ 2 ■ s^ 1 /K,AT 1 ~ 

10~ 4o K and H = 10cm we obtain for A$ G (Gauss ■ cm 2 ) ~ 3 x 10~ 19 • Ar(s). Measuring 

such tiny variations of magnetic field for large At is near to the limit of SQUID sensitivity. 

The main problem in observation in the flux (21) and (22) will be connected with gener- 
ating thermal current J s = —A~ 1 / 2 \f3gradT comparable in magnitude with its critical value 
j c , since the temperature difference across the junction is limited by a low temperature T c 
and small sizes of the junction. 

The similar method of measurement has been used by Jain et al [4] in null result exper- 
iment on confirmation of the strong equivalence principle for a charged massive particle. In 
their experiment the phase of Josephson contacts has been locked to an external microwave 
source and is schown to be technically feasible to measure a voltage drop 10~ 22 V. 

It is interesting to mention that the predicted mechanism for production of magnetic 
field and current changing with time can be of crucial importance in astrophysics as a way 
(additional to the existed ones [17]) for generation of electromagnetic radiation from pulsars. 
According to the recent theoretical models [7], a neutron star is the relativistic compact 
object consisting of the conducting crust and superfluid core. In the inner crust of the 
neutron star the superfluid coexists with a crystal lattice and in its core, at densities above 
2 x 10 1A gm/cm 3 there is a homogeneous mixture of superfluid neutrons and superconducting 
protons. 

Important fact is that the thermoelectric power (3 is the function of temperature as T 3 / 2 
and in this connection can reach large numbers since superconductivity in the stars takes 
place at the temperatures 10 6 -i-10 7o -ft'. So if one accepts that the SNS structures are realized 
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in the intermediate boundary between conducting crust and superconducting core inside the 
neutron star then the strong heat fluxes in these SNS junctions can lead to the generation of 
time-varying magnetic field (i.e. electromagnetic radiation) due to the thermoelectric effect 
described by the basic formula (21). 

Suppose that, in addition to the heat flow, SNS Josephson junction is placed in a 
magnetic field parallel to the plane of the junction (the xy plane). It is well-known that the 
maximal current density, j( s ) ma xa, which can pass through the junction is 

_ A sin7r£b/£o , 9 „v 

J(s)maxa J(c)a\ 7r <$> b /<$> Q ^ ' 

where $6 is the magnetic flux through the junction. 

Taking into account that j( s ) a = A/3A -1 / 2 9 Q T we can find that the critical temperature 
difference, corresponding to appearance of voltage across the junction, is 

(Af) c = A(0)^|!^^| (24) 

and depends from gravitational field, where A = A(0) at the junction. This is the general 
relativistic generalization of a thermal analog of the dc Josephson effect [18,19] according to 
which the critical value of heat flow through a SNS junction is a nonmonotonic function of 
the magnetic flux <3> b . 

Thus the gravitothermoelectric phenomena in superconductors considered here allow us, 
in principle, to detect the general relativistic effects. But we would like to emphasize that we 
only concentrated on gravitothermoelectric phenomena in superconductors of type I. Never- 
theless, recently several thermoelectric effects were observed in high temperature supercon- 
ducting materials of type II (see, for review, [5]). In this connection further investigation is 
needed for taking into account gravitational corrections for thermoelectric effects in type II 
superconductors. 
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